Small parameter for lattice models with strong interaction 
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Diagram series expansion for lattice models with a localized nonlinearity can be renormalized so 
that diagram vertexes become irreducible vertex parts of certain impurity model. Thus renormalized 
series converges well in the very opposite cases of tight and weak binding and pretends to describe 
in a regular way strong-correlated systems with localized interaction. Benchmark results for the 
classical O(N) models on a cubic lattice are presented. 



One of the key problems of the modern condensed- 
matter theoretical physics is the quest for a regular de- 
scription of the strong correlations [l| . Roughly speaking, 
correlations are strong if the nonlinearity and coupling 
are comparable, and therefore perturbative approaches 
make no sense. 

Consider a lattice model with local coupling between 
cells with a nonlinearity localized inside each cell. In the 
two limit cases of the weak and tight binding a series 
expansion can be performed in powers of the respective 
small parameters |2| ■ One of a few regular ways to han- 
dle the situation of crossover is to employ the concept of 
localized correlations. It basically means that interaction 
of each cell with the rest of the system can approximately 
be described as an interaction with certain non-correlated 
effective bath The parameters of this bath are to 
be found self-consistently. If there is, say, one particle 
per cell, it means that real multi-site problem is approx- 
imated by much simpler single-site one ("impurity prob- 
lem" ) . An example of such an approach is the mean-field 
description of the classical statistical ensembles of inter- 
acting particles. In this scheme, interaction of a particle 
with its fluctuating surrounding is replaced with inter- 
action of a particle with a constant average field. There 
is a self-consistent equation for this field. An accuracy 
can be dramatically increased for change the mean-field 
by Gaussian-bath approximation. In this case, the sur- 
rounding is approximated by certain Gaussian ensemble. 
It is important that almost the same argumentation can 
be presented for very wide class of systems. Particu- 
larly, so-called dynamical mean-field theory (DMFT) for 
strongly interacting fermions is essentially a Gaussian- 
bath approximation 

Physically, Gaussian-bath approximation consists in 
the assumption that the irreducible part of high momenta 
is localized inside a single site. The strong point of this 
assumption is that it is valid for several very different lim- 
its for systems with localized nonlinearity. The approxi- 
mation reproduces leading-order results of the weak- and 
tight-binding series expansion simultaneously. Indeed, if 
the nonlinearity is small, the perturbation theory explic- 
itly says that correlations are localized in the first-order 
of weak-binding approximation. On the other hand, if 
coupling goes to zero, sites become almost isolated and 
correlations are obviously localized, no matter how strong 



the nonlinearity is. Once the two opposite limits are de- 
scribed well, their crossover can be expected to be also 
somehow depicted. There is some more justifications for 
the scheme, particularly it becomes exact for a system of 
an infinite dimensionality. 

Since an analysis of proper impurity problem gives 
a good understanding of the properties of strongly- 
correlated systems, it is desirable to construct a pertur- 
bation theory starting from the solution of an impurity 
problem as a zeroth approximation. A theory of this 
kind is presented in this paper for lattice models with a 
localized nonlinearity. We exactly renormalize diagram 
expansion so that diagram vertexes become irreducible 
vertex parts of certain impurity problem. An impor- 
tant peculiarity of the proposed approach is that a good 
convergence is achieved both in tight-binding and weak- 
binding situations. As an example, we present the results 
obtained for quantum Ising model in a transverse field. 

We consider statistics of the real vector field on a 
discrete lattice with a localized nonlinearity 0(0) and 
the dispersion law f2. In general, lattice potential can 
depend on a number of the lattice site. Through the 
article we number sites by x whereas components of the 
on-site field 4> x are numbered by j, statistical average is 
denoted by triangle brackets: 

<...>= Z- 1 J ... e-^^iDcj)} (1) 

Z = f e -f>EW[ D Q 

^) = Ex9x(W + i(W). 

Scalar product (</>Sl(/>) is a sum over both indices: 
Tlixx'jj' 4 > xjQxx'jj' ( t > x'j'] below we also use external prod- 
uct denoted by the dot: <j) ■ cf>. The notation {4> x Qxx4>x) 
is used for a scalar product in the on-site subspace: 

{^(j^x^xx^x) = ^xj^xxjj'^xj'- 

Consider an auxiliary ensemble with energy 

E avx {4>) = ]T Qxi^x) + \ (<M0) + (4,°(n -A)4>). (2) 

X ~ 

We introduced here vector quantity (f>° and a Hermi- 
tian tensor A having non-zero terms only at the same 
x-indices: A xx ' = for x =/= x' (through the text we call 
such quantities 'ir-diagonal'). The values of A and <p° are 
not specified for a while. 
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Since A is x-diagonal, site oscillators of the auxiliary 
ensemble are uncoupled (one can see that E aux ((f>) splits 
to a sum of independent single-site energies). Thus, prop- 
erties of the auxiliary system at each site can be obtained 
as a solution of a single-site impurity problem. 



Denote an average over the auxiliary ensemble by over- 
line and introduce normalized two-point correlator and 
higher-order irreducible vertex parts at the sites of aux- 
iliary system: 



.912 = (3 (4>i4>2 - 4>i 02 ) (3) 

7i?2'3' = -P 3 9l'l92>293>l (010203 ~ 0102 03 ~ 0103 02 ~ 0203 01 + 201 <t>2 fo) 



7i' n> = ( — 0) n 9vi—9n'n ' (irreducible part of the n— th moment). 



To simplify the notation, we omit site index x and write 
numbers as subscripts for the on-site coordinates here. 

It should be noted that mean-field and Gaussian-bath 
approximations can be understood as approximations of 
the system (QJ by (J3J) with a self-consistent requirement 
— (f> and certain choice of A xx . Indeed, to obtain the 
mean-field scheme, one should, first, replace all <f> in E{<j>) 
with its average for all cells except just a single one, and, 
second, use an average over this particular cell as a guess 
for < <f> >. It is easy to check that the mean- field corre- 
sponds to the simple choice A xx — Q xx . Gaussian-bath 
approximation also reduces the system to a single-site 
problem. The difference is that the remaining degrees 
of freedom are not frozen in its average positions, but 
approximated with harmonic oscillators, forming Gaus- 
sian bath. These Gaussian degrees of freedom can be 
integrated out, giving again single-site energy like in |j2J 
with certain A. An approximation for the single-site sys- 
tem by a harmonic oscillator can be found from the re- 



quirement that the approximation mimics <f> x and g xx . 
Thus obtained parameters of harmonic oscillators are to 
be used as in the next iteration of the self-consistent loop 
as parameters of the oscillators forming the bath. These 
iterations converge to the point where s-diagonal part of 
the tensor (g + (fi — ^4) _1 ) vanish. This is in fact a 
condition for A of the Gaussian-bath approximation. 

The aim of this paper is to express averages Q via 0> 9 
and 7'"' by certain regular series. The quantities A and 
(po will determined by certain self-consistent conditions. 

First, we would like to introduce the "dual" variables /. 
Let us consider the basis where g(fl — A)g is diagonal and 
label the states in this basis by p. Define the integration 
J[Df] = J df Pl ... J df Pn ... in such a way that the path of 
integration in each integral depends on a sign of diagonal 
element (g(fl — A)g) pp : f p goes from —00 to 00 for (g(Cl — 
A)g) pp < and from —ioo to ioo overwise. With thus 
defined J[Df], the following identity holds: 



±p{(<p-<p a )(n-A){<p-<p a )) = ( factor ) x f [ jD j] e -/3(((0-0°) 3 - 1 /)-|(/(5(^-^)s)" 1 /)), (4) 



In fact, the path of integration over f p is chosen to deliver the convergence of integrals here. After this identity is 
applied, partition function takes the form 

Z = f[D4>]f[Df\e-^f) (5) 
E(<P, f) = E—{4>) + {{4> - O )<T 1 /) - |(/(ff(n - A)g)- 1 f), 

A set of exact relations between the averages over initial and dual ensemble can be established by the integration 
by parts in JSJl with respect to /. Particularly, 

g(Q — A) <4> — 4>° >=< f >; (6) 
g(Q - A) < (4> - cfP) ■ (0 - 0°) > (n - A)g =</•/> - A)g. 

I 



The idea of an introduction of new variables is that energy (JSJ) does not contain direct coupling between dif- 
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FIG. 1: Diagrams contributing E' and corresponding coeffi- 
cients. Diagram (a) is small if condition (11) is fulfilled. Di- 
agrams (b) and (c) describe av a 2 correction to E'; diagram 
(c) vanishes for an unordered state of O(N) model. 

ferent sites for the initial variable 0, since (</>3 -1 /) — 
Ei(^(9 _1 /)i;)- Therefore <j> x can be integrated out at 
each site, yielding an an expression for the energy in dual 
variables E(f): 

E(f) = EMu) + k(m) (7) 

n = -g- 1 -g- 1 (n-A)- 1 g- 1 

This expression is formally very similar to an initial sys- 
tem QJ. Physical difference comes, first, from the pos- 
sibility to choose A and <f>° in certain optimal way, and, 
second, from the the non-trivial way of integration over / 
so that / has an imaginary part. In particular, < /J - > 
can be negative or equal zero. 

Properties of the dual potential V(f) are determined 
by the cumulant expansion for an auxiliary system. 
The last term in JJJ is chosen in such a way that the 
second derivative of V vanishes; other derivatives are 
= - (jP)g- 1 and -pV™ = 7 (n) for n > 3. To 
describe effects due to nonlinearity, we introduce tensor 
E' defined by equation 

<0.0>= p- x {g- 1 +n- A + E')~\ (8) 

so that £' describes a non-local contribution to the self- 
energy of the initial system . It follows from the second 
line of © that dual two-point correlator can be expressed 
via £ as follows: 

< f ■ f >= %-i-v)~-?+(n-A)-i 9 - (9) 

All expressions presented above are exact. Now, we 
construct a perturbation theory resulting from the series 
expansion in powers of V. Terms of these series can be 
expressed by the diagrams with standing at vertexes 
connected with lines carrying dual two-point correlator < 
/•/ >. We consider the series for E'. It can be found from 
the expansion of the left- and right-hand side of © in 
powers of V and E', respectively. Some of those diagrams 
are presented in Figure 1. Each diagram is accompanied 
by an additional factor — P^ 1 and a numerical coefficient. 
It can be found from the above-mentioned expansion of 

©■ 



Formally, single-lag vertexes can occur in the diagrams 
due to the presence of in the dual energy J7J. We 
require such (f>° that < / >= to get rid of them. Then 
the first line of JJyl requires 

< >= 0°. (10) 

Since 7^™) of any order can in principle appear in the 
series, a formal small parameter should be introduced 
to group the terms properly. Consider a special case 
of a scalar field with small nonlinearity: U{4>) — a0 4 , 
a — > +0. Estimate how a appears in high derivatives 
of V(f): 7( 2 ™) oc a n and -y^ 2 ™ -1 ' oc a n <f>o- We propose 
to group diagrams with respect to their formal small- 
ness in powers of a. This is clearly a good choice for a 
weak-binding limit. Consider the opposite case of tight- 
binding £7 — ^ and suppose that our choice of A delivers 
also (A — fi) — > (note that mean-field and Gaussian- 
bath schemes satisfy this condition). It is obvious from 
that bare two-point correlator appears to be a small 
quantity in this case. Since the number of lines in a 
diagram is roughly proportional to its a-order, a again 
appears to be a good small parameter. We conclude that 
series expansion converges well both in tight-binding and 
weak-binding limits and therefore expect that the theory 
is suitable for a crossover situation. 

Formally, series expansion can be constructed for any 
A, but an appropriate choice should result to better con- 
vergence. It is useful to require 

r x g =< (0 - 0o) • - <h) >xx ■ (11) 

One can formally obtain Gaussian-bath approximation 
as a result of Gaussian approximation for < / • / > in 
©, combined with requirements i|l(Jllll) . Gaussian ap- 
proximation means a order of the theory. Higher terms 
of an expansion in powers of a improve the result and, 
in particular, bring non-local correlations on the scene. 

To demonstrate capabilities of the method and to give 
more technical details, we present here the results ob- 
tained in an a 2 -approximation for classical O(N) models 
at 3D cubic lattice with the nearest-neighbor coupling. 
Potential energy is 

<xx'> 

where <j) x is an TV-component vector of the constrained 
length I \<f) x 1 1 2 = (4> x (f>x) = N. In the previous notation, 
dispersion law is flk — ~ (cos fc^+cos fc^+cos k z ), whereas 
the constrain makes the system nonlinear. Sum is over 
all pairs of the nearest neighbors. Cases of N = 1,2,3 
correspond to Ising, xy, and Heisenberg models, respec- 
tively. In three dimensions, system (|12|l shows a second- 
order phase transition at finite temperature for N. The 
model is extensively studied. One can note Monte Carlo 
simulations 0, high-temperature expansion (in 
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our terminology, this is in fact tight-binding approach), 
renormalisation-group methods [8j etc. Here we do not 
pretend to obtain something new for the model itself. It 
is used to check if the method works for a realistic system. 

For simplicity, let us consider an unordered state 
< <j) >= 0, so that only even-lag vertexes appear 
in diagrams, 'y( 2n + 1 ) =0 . It can be easily obtained 
that 512 = 0~ l &i2 (irrespective on N) and 7^34 = 

77v(<5l2<534 + ^13^24 + ^14^23) with 7jv=1 = -2/3, 7jv=2 = 

— 1/2,7n = 3 = —2/5. Here <5 is a Kronecker delta. Ten- 
sor structure of g, 7^ 4 \ as well as A, corresponds to the 
O(N) symmetry of the problem. 

The self-consistent equation for a follows from |(HJ and 
CD: 

If d 3 k 
9 ~ (2^)3 J 5 -i + il k + £' - A ' ' 

Technically, this equation can be solved iteratively, 
similarly to DMFT-loop technique p|. Start from cer- 
tain guess for A and £'. For a general case, one should 
calculate g(A) and j( n >(A) at this point, but for a par- 
ticular model i|12|) it is not necessary as g and 7 do not 
depend on A, because of the constrain ||</>x|| 2 = N. Cal- 
culate < f ■ f > from (J5J, given g, A, and £'. For thus 
obtained < /•/ > and j( n ', calculate new guess for £' for 
a given order in a (see below expressions for £' up to a 2 ). 
Finally calculate the right-hand side of 1)13(1 . given g, A, 
and new Denote the obtained value by g and take 
the quantity A + ^g^ 1 — g~ l ) as new guess for A, where 
£ < 1 is a numerical factor. Repeat the loop for new A 
and £'. For properly chosen initial guess and value of 
£ iterations converge to a fixed point satisfying equation 

Gaussian-bath approximation corresponds to an as- 
sumption £' = 0. Corrections are given by the dia- 
grams presented in Figure 1. Lines in these diagrams 
are thick, i.e. they correspond to the renormalized cor- 
relator < / ■ / >. First-order correction would be given 
by a simple loop shown as diagram (a) in Figure 1. This 
diagram corresponds to formula — ^g7^ 4 - ) < fx ' fx >■ 
It is however absent, because < f x ■ f x > equals zero in 
Gaussian-bath approximation, as it can be obtained from 
© and (|13|l . So the correction starts actually from a 2 
terms. 

Strictly speaking, once high-order corrections are taken 
into account, < f x ■ f x > does not vanish anymore. But 
it remains small. It can be shown that simple loops does 
not appear in diagrams up to a 3 order. Consequently, 
second-order correction to £' can be described by just 
a single diagram (b) drawn in Figure 1, so that £^. x / — 
~ ftp^ 2 ^ f x ' ^ x> >3 ' After the sum over internal indices 
is taken, it gives T,' xx , = -^"VU + N/2)fl x ,, where 
f x ,x' is a diagonal clement of the dual correlator. 

Figure 2 shows data for the effective soft-mode dis- 
persion £0 = (P < 4>k=o ■ 4>k=o >) • Dots are numer- 
ical data. Data for the lines are obtained by formula 



So = g^ 1 + ^fc=o + £' — A. Iteration procedure converges 
well everywhere except very narrow regions near critical 
point. Since we were not interested in a description of 
the critical point itself, we did not put special attention 
to the source of this divergence. 

Gaussian-bath approximation predicts the same de- 
pendence for any N, as one can observe from ((13(1 with 
£' = 0. This is a serious qualitative drawback, although 
the quantitative accuracy is not so bad. So, the depen- 
dence of £0 on N in the model is related with a non-local 
part of correlations, and should be described by higher 
terms of a-corrections. Indeed, a 2 -correction describes 
the dependence T,q(N) well. There is also a dramatic 
increase of the accuracy in the a 2 curves. 

For the broken-symmetry state, the theory can be con- 
structed in a very similar manner. Iteration loop now 
includes also a modification of <jp to fulfill the condition 
(|1U[1 at a fixed point. Diagram (c) from Figure 1 appears 
due to the presence of 7 < - 3 - ) . We do not show numerical 
results here, because they are very similar to what we 
obtain for the unordered state. 

The entire plot range of Figure 2 lies within a critical 
region. For example, all points for Ising model obey crit- 
ical scaling law < ^ 2 . =0 >cx (/3 — /3 C ) 124 with a good accu- 
racy, as the dashed line shows. Thus the theory performs 
well deep inside the critical region. However, the very 
vicinity of the critical point is not correctly described, 
because we work with perturbation theory of a finite or- 
der. It would be worth to construct a renormalization- 
group theory based on the presented perturbation ap- 
proach. An attempt to construct a primitive approach of 
this kind starting from mean-field theory was performed 
several years ago [loj. 

There is a comment about an applicability of other 
approximation schemes. The system is strongly nonlin- 
ear, so weak-binding approximation is clearly not valid 
here. On the other hand, tight-binding description of a 
low order fails near the phase transition point because 
of the critical increase of the correlation length. For 
example, second order of the tight-binding series gives 
£ = /3 _1 (l+4/3) _1 . This is very inaccurate; correspond- 
ing curve simply lies out of the plot range of Figure 2. 
More sophisticated schemes, for example, cumulant ex- 
pansion, also require much higher order to achieve an 
accuracy comparable with the presented a 2 results. It 
is also important to recall that Gaussian-bath approxi- 
mation becomes exact in the limit N — > 00, and X/N 
expansion can be constructed at this point 11]. But 
1/N approach is essentially based on a very particular 
O(N) symmetry of system 1(12(1 . Contrary, our method 
is developed for a general case and does not require any 
special symmetry. 

In conclusion, we presented an exact renormalization 
of the diagram-series expansion in terms of the self- 
consistent impurity model. Renormalized vertexes are 
irreducible vertex parts of the impurity model. There 
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FIG. 2: Effective soft-mode dispersion for O(N) model at TV = 1,2,3. Dashed line: Gaussian-bath approximation. Solid lines: 
results for the a 2 -scheme, from bottom to top: N — 1 (Ising), N = 2 (xy), and N = 3 (Heisenberg). In a narrow region near 
critical point the iteration procedure did not converge, and the lines are not drawn. Symbols are numerical data for TV = 1, 2, 3: 
respectively solid circles, open circles and crosses. Dotted line is a fit corresponding to a 'purely critical' behavior of < <f) ■ <j> > 
in Ising model. 



is an explicit small parameter in the theory for weak- 
binding and tight-binding limit. A worked example of 
an a 2 -approximation for simple 3D models demonstrates 
that the method performs well for a strong-correlated sit- 
uation. The example models were chosen because of the 
physical simplicity; the method itself looks rather gen- 
eral. It would be particularly important to extend it 
to the systems of identical quantum particles. Another 
challenge is a construction of renormalization-group the- 
ory based on the method. 
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